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ABSTRACT 

We study boundary reflection matrix for the quantum fleld theory deflned on 
a half hne using Feynman's perturbation theory. The boundary reflection matrix 
can be extracted directly from the two-point correlation function. This enables us 
to determine the boundary reflection matrix for affine Toda fleld theory with the 
Neumann boundary condition modulo 'a mysterious factor half. 
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I. Introduction 



The S'-matrix of (massive) integrable quantum field theory in 1+1 dimensions can 
be studied by several different methods. The high brow technology to construct S'- 
matrix is based on the symmetry principle such as Yang-Baxter equation, unitarity, 
crossing relation, real analyticity and bootstrap equation|I|, 0, This program 
entirely relies on the assumed quantum integrability of the model and produces an 
S'-matrix which is exact up to all loop order. 

Despite its beautiful nature, this method inevitably needs additional information. 
Furthermore, there is an inherent so-called CDD ambiguity. To back up the situation, 
Feynman's perturbation theory has been used and shown to agree well with the 
conjectured 'minimal' S'-matrices[0, ||, H, ||, |^, ||, P]. This may also be considered 
as the strong evidence for the assumed quantum integrability. In perturbation theory, 
S'-matrix is extracted from the four-point correlation function with LSZ reduction 
formalism. 

About a decade ago, integrable quantum field theory on a half line (—00 < x < 0) 
was studied using symmetry principles under the assumption that the integrability 
of the model remains intact |jlO[. The boundary Yang-Baxter equation, unitarity 
relation for boundary reflection matrix K^{0) which is conceived to describe the 
scattering process off a wall was introduced [PH]]. 
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Figure 1. Boundary Reflection Matrix. 

Recently, boundary crossing relation was introduced [|ll|]. In fact, the bound- 
ary crossing relation is automatically satisfied if the boundary bootstrap equation 
Subsequently, some exact boundary reflection matrices have been 



is satisfied [12 



constructed ITT], |T3|, |T^, |T4|, [T^ for afiine Toda field theory(ATFT). However, it turns 
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out that there is a plethora of solutions for boundary reflection matrix despite re- 
quiring the 'minimality' assumption which has been effective in the S-matrix theory 
on a full line. Furthermore, it is unknown what particular boundary condition(or 
boundary potential) actually corresponds to a particular solution. 

In order to have a direct access to the boundary reflection matrix, it seems 
compelling to study the boundary system from the Lagrangian quantum field theory. 
In fact, several studies on the boundary system has already been done|26|, 0, |l^ 
in the Lagrangian quantum field theory context. However, the boundary reflection 
matrix has never been discussed in this framework, yet. 

On the other hand, the ordinary LSZ reduction method which does the job to 
extract the S-matrix from the off-shell correlation functions becomes inapplicable 
for the quantum field theory on a half line, since the momentum eigenstates with a 
definite sign in the asymptotic region do not satisfy the boundary condition at the 
origin in space. 

In this paper, we propose a method to extract boundary reflection matrix directly 
from the two-point correlation function. In section II, we describe the formalism. 
In section III, we present the one loop result for the sinh-Gordon model (or aj^^^ 
affine Toda theory). In section IV, we present the one loop result for the Bullough- 
Dodd model (or 02 affine Toda theory). We also give a conjecture for the exact 
boundary reflection matrix guided from this one loop result. This model fully utilises 
all possible Feynman diagrams because it has three point self-coupling as well as 
four-point coupling. Finally, we give some discussions in section V. 

II. Boundary Reflection Matrix 

We are mainly concerned with affine Toda field theory with integrable boundary 
interaction though the formalism may be applicable for any quantum field theory. 

To begin with, we review the two-point function for the model on a full line. 
The bosonic ATFT]^ is defined by the following Lagrangian density based on a Lie 
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algebra g with rank r. 

£($) = \d,rd^<t>^ -"if: ^^e'^""*' (1) 

where 

r 

ao — — ^njOJi, and no = 1. 
1=1 

The field 0" (a = l,---,r) is a-th component of the scalar field and ctj {i — 
1, • • • , r) are simple roots of g normalized so that the universal function B{(3) through 
which the dimensionless coupling constant (3 appears in the (S'-matrix takes the fol- 
lowing form: 

- ijTTWWy 

The m sets the mass scale and the n^s are the so-called Kac labels which are char- 
acteristic integers defined for each Lie algebra. 

The two-point function at tree level is given by the Feynman propagator: 

G{t', x'- t,x)^ [ .% , e-(^--^)+-^(-'-), (3) 

J [277)^ — + le 

where p = (w, k) is the two dimensional energy-momentum and is the mass of the 
particle in the original Lagrangian. As is well known, this two-point function depends 
only on the difference of its arguments and accommodates contributions coming from 
the positive energy states as well as the negative energy states depending on the sign 
of the difference of the time arguements. This physical feature is usually implemented 
by is prescription or choice of w-contour. 

At one loop order, there are three types of Feynman diagram contributing to the 
two-point correlation function as depicted in Figure 2. 



4 





Type I. 



Type II. 



Type III. 



Figure 2. Diagrams for the one loop two-point function. 

Type I, II diagrams have logarithmic infinity independent of the external energy- 
momenta and are the only divergent diagrams in 1+1 dimensions. This infinity is 
usually absorbed into the infinite mass renormalization. Type III diagrams have 
finite corrections depending on the external energy-momenta and produces a double 
pole to the two-point correlation function. 

The remedy for these double poles is to introduce a counter term to the original 
Lagrangian to cancel this term(or to renormalize the mass). In addition, to maintain 
the residue of the pole, we have to introduce wave function renormalization. Then 
the renormalized two-point correlation function remains the same as the tree level 
one with renormalized mass m^, whose ratios are the same as the classical value. 
This mass renormalization procedure can be generalized to arbitrary order of loops. 

Now we consider the model on a half line(— oo < x < 0). The action is defined 
as follows. 



/O foo ( 1 jTn2 r \ 

-oo J—oo jj y 



(4) 



The above simple action may have additional boundary potential maintaining the 
integrability. Non-trivial boundary potentials which do not destroy the integrability 
have been determined at the classical level [ITTl, [ 
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l[|. The stability of the 
. Here we consider the 



model with boundary potential has also been discussed [|n 
model with no boundary potential, which corresponds to the Neumann boundary 
condition: ^7 = at x = 0. This case is believed to be quantum stable in the sense 
that the existence of a boundary does not change the structure of the spectrum. 
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At tree level, two-point correlation functions are given by a sum of a direct 
contribution and a reflected one which may be considered as coming from the image 
point, 

GN{t',x']t,x) = G(t' ,x';t,x) + G(t' , x';t, —x) (5) 



j2 

P ^ ^~iw(t'~t)(^ik(x'-x) j_ ^ikix'+x)-^ 

(2'nY — + is 
We may use the /c-integrated version. 

GM{t\x';t,x) = f ^}^^-Mt'-t)^^r-kW-x\ ^ ^-fkix'+x)^^ ^ I 2 _ ^2. (Q) 

J 111 Ik, * 
We find that the unintegrated version is very useful to extract the asymptotic part 
of the two-point correlation function far away from the boundary. 

To compute two-point correlation functions at one loop order, we follow the 
idea of the conventional perturbation theory p!6|, |T^. That is, we generate the 
relevant Feynman diagrams and then evaluate each of them by using the zero-th 
order two-point function for each line occurring in the Feynman diagrams. 

Type I diagram gives the following contribution. 

roo 

dXi / dtiG iy(t, X] ti, Xi) Gjyit ,X ; ti, Xi) Giy{ti, Xi] ti, Xi). (7) 



oo 



Let us take a close look at GN(ti,xi;ti,Xi). 

= f^- ^^--(l + e^'^^^-). (8) 

J {2iT)^ pf — mf + le 

The first term is the ordinary infinite mass renormalization term as for the full line 
theory. To cancel this, we introduce a counter term exactly the same as for the 
full line. We should not simply discard the second term and this term contributes 
to the boundary reflection matrix. We evaluate ti integral in Eq.(|^, giving energy 
conservation at the interaction vertex. The Xi integral gives 'spatial momentum 
conservation' as follows: 

k ±k' + 2ki = 0. (9) 

After the integrations of loop variables and energy conserving delta function 
which resulted from ti integral, we get the following result from Type I diagram, 

r dw^dk^dl^ g-ii«(t'-t)gi(fca;+fc'x') i i J QgN 

J 2tt 2n 2tc w"^ — k"^ — + ie — k'"^ — m\ + ie ' 
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where ki is defined in Eq.(^. 

From Type II diagram, we can read off the following expression: 

rO 



dxidx2 dtidt2GNit, x;ti, xi) G^it' , x';ti, xi) GNih, Xi;t2, X2) (11) 

D J — OO 

G]\f(t2, X2; t2, X2). 

Similarly as for the Type I diagram, GN{t2,X2',t2,X2) contains the ordinary infinite 
tadpole term. By introducing infinite mass renormalization, we can discard this 
tadpole term. 

The ti and ^2 integral give energy conservations at each vertex. The xi and X2 
integral gives 'spatial momentum conservation' as follows, 

k ±k' + ki = 0, ki + 2k2 = 0. (12) 

After the integrations of loop variable W2, momentum conserving delta functions 
as in Eq.([I^) and energy conserving delta functions 6{wi),6{w' + w), we get the 
following result from Type II diagram, 

dw^Md^ ^_iw(t'-t)^i{kx+k'x') i i j^23) 

271 271 271 w"^ — k'^ — ml + iew"^ — k'"^ — ml + ie ' 



H + mi 2^kl + ml 
Type III diagram gives following contribution: 



dxidx2 dtidt2GN{t,x;ti,Xi) G^it' ,x';t2,X2) GN{t2,X2;ti,xi) (14) 

J—oo 

GN^h, X2; ti, Xi). 

The ti and t2 integral give energy conservations 6{w + wi + W2), S{wi + W2 — w') 
at each vertex. The xi and X2 integral give 'spatial momentum conservations' as 
follows, 

k ±ki±k2 = 0, ±ki ±k2 + k' = 0. (15) 



Among the 16 possible combinations of signs in front of each spatial momentum in 
the above equation, the same 8 combinations as the momenta conservation on a full 
line gives exactly the same finite mass renormalization. The other combinations give 
the following result from Type III diagram, 



/ 



dwdkdk^ ^-iyjit'-t)^i(kx+k'x') i i /^gx 

27r 27r 27r w'^ — — ml + iew'^ — k'^ — m^ + is ' 

I ~^ { ' + ' ) 

^ 4 2wi{wi — Wi){wi — Wi) {wi — Wi){W]^ + Wi){W]^ — Wi) ' 

where we included | which was introduced while we were extending the domain of 

xi,X2 integrations to a full line to allow the delta function and we introduced the 

following notations. 



\Jkl + ml, wf^w + ^Jk2 + ml, w-^ ^ w — \fl^^^m!i. (17) 



It should be remarked that this term should be symmetrized with respect to m^, mc- 
Now we propose a method to extract boundary reflection matrix directly from 
the two-point correlation function. The general form of each contributions coming 
from type 1,11 and III diagrams can be written as follows: 

d^dkdk^^-Mt'-t)^iiM'.') J J 7(^, k, k'). 

2tt 2n 2tt — k^ — m^ + ie — k'^ — m^ + is 

(18) 

Contrary to the other terms which resemble those of a full line, this integral has two 
spatial momentum integration. First, let us consider the k' integration. There are 
two contributions. One comes from the usual pole contribution of the propagator 
and the other one from the poles and the branch cuts of / function if any. For the 
k integration, the similar consideration can be done. Here we simply state that the 
contributions other than the usual pole contributions coming from each poles of the 
external propagators turn out to be exponentially damped as x,x' go to —oo. 

In this way, we can get a method to compute elastic boundary reflection matrix 
Ka{0) deflned as the coefficient of the reflected term of the exact two-point correlation 
function in the asymptotic region far away from the boundary. 

f 1 ^ i^^(^)e-S(-'+-)), k = ^w^ - ml 

J ATT Zitv 
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Ka{6) is obtained using w = niaCoshO. 

Here we list each one loop contribution to Ka{6) from the three types of diagram 
depicted in Figure 2: 

^i'H^) = ttAtttI ; . i„ +7^)CiS„ (19) 



2mashe 2jmWe + m? 2m, 



6 

— Z 
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= 7^ ^( ^ , + TT^) ^2 ^2, (20) 

^^'''H^) = 7^^(4/3(A;i = 0,fc2 = A;)+4/3(A;i = fc,A;2 = 0))C3^3, (21) 

ZTTlaSllU 



where I3, is the function defined in Eg. ([Tq) , and the factor 4 in front of it accounts 
for the fact that there are four combinations in Eq.(^) which give the same result. 
The Ci, Si denote numerical coupling factors and symmetry factors, respectively. 

III. Example I : ai^^ affine Toda theory 

For the sinh-Gordon model, only Type I diagram in Figure 2 contributes to one loop 
two-point correlation function since there is no three-point self coupling. We have 
to fix the normalization of roots so that the standard -B(/3) function takes the form 
given in Eq.(|D. 

We use the Lagrangian density given as follows. 

1 

— ( 

2 



£(0) = -9^09'^0-V(0), (22) 



4/32 



^(0) = ^(eV^/3^ + e-v^'3^-2) (23) 



^mV-^mW + 0(/3^) 



The scattering matrix for the elementary scalar of this model is|l|] 

Here B is the same function defined in Eq.(^ and we used the usual notation of 
building block[0] as follows. 

_ sh{e/2 + z7rx/2/^) 
^' sh{e/2-tnx/2h)- ^ ' 
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For the sinh-Gordon model, h = 2 and from now on we set m = 1. 
The result coming from Type I diagram is 

It turns out that this is too large by a factor 2 to satisfy the crossing unitarity 
relation at one loop order. 

K{e) K{e -tn) = s{2e). (27) 

So we need to include an extra factor | into our formulae in Eg . (|2l|) , although 
we do not understand the reason. Then, we find that the formulae in Eq. (|2lD with 
the extra factor ^ work for any theory. 

On the other hand, there are two 'minimal' boundary reflection matrices which 
are meromorphic in terms of rapidity variable are known for ai'^ modelO, KM. One 
of them agrees with the perturbative result. 



where 



K{9) = [1/2], (2^ 



fd = (x-l/2)(x+l/2) 

{x-l/2 + B/2){x + l/2-B/2)' ^ ' 



(2) 

IV. Example II : a2 affine Toda theory 

The BuUough-Dodd model has three-point self coupling as well as four-point cou- 
pling. So all possible three types of diagram in Figure 2 contribute to one loop 
two-point correlation function. We have to fix the normalization of roots so that the 
standard function takes the form given in Eq.(0). 
We use the Lagrangian density given as follows. 



1 

2' 

2 

m 



C-i'P) = ^d^<pd^<p-V{4>), (30) 



V{ct>) = -^(2e^<^ + e-2^'^-3) (31) 

1 0,0 1 o-^,*? 1 



' 6 8 
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The scattering matrix for the elementary scalar of this model isQ 

^^^>~ {B){2-B){l + B){3-By 

For the BuUough-Dodd model, h = 3. 

The result coming from Type I diagram is 

The result coming from Type II diagram is 

For type III diagram, when ki = 0,k2 = k, 

= 1, w+ = 2che, = 0, (35) 

and when ki = k,k2 = 0, 

Wi = ch9, = ch9 + 1, = ch9 — 1. (36) 

The result coming from Type III diagram is 

j^(ni) = J_(^ 1 + i + _J_ + 1 ) 

ishe^ 2{1 - 2ch9) {2che - l){2chd + l)2che -2che {2che + 1)2^ 

(37) 

x(^/3)2xl8. 



Adding above three contributions as well as the tree result 1, we get 

, shO shO shO sh6 shO rt( f^^\ 

K{0) = + ^(-^^^^ - 2c/,^ - V3 ^d^^ che + 1/2 " 2che + ^3^ ^ ^' 

(38) 

This satisfies boundary crossing unitarity relation and boundary bootstrap equation 
up to order. 

K{e) K{e -in) = s{2e), K{e) = K{e + iTi/3) K{e-m/3) s{2e). (39) 

On the other hand, there are four 'minimal' boundary reflection matrices which 
are meromorphic in terms of rapidity variable are known for 02 model [14, 23 1. None 
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of these corresponds to the perturbative result. A possible exact solution would be 
the following. 



This is one of the 'minimal' non-meromorphic solutions with square root branches 
which are determined from symmetry principles such as boundary crossing relation 
and boundary bootstrap equation. 

V. Discussions 

In this paper, we proposed a method to compute boundary reflection matrix directly 
from the two-point correlation function rather than using the LSZ reduction which 
is not applicable to the quantum field theory on a half line. In our formalism, the 
unintegrated version of the Neumann Green's function turns out to be very useful to 
extract the asymptotics of the two-point correlation function in the region far away 
from the boundary. This enables us to determine the boundary reflection matrix 
for the affine Toda field theory, specifically o^l^ and 0,^2^ models, with the Neumann 
boundary condition modulo 'a mysterious factor half. 

We have also done similar computations for some a^^^ models as well as some d^^^ 
models. When the theory has a particle spectrum with more than one mass, each 
contribution from three types of diagram in Eq. (|2l|) has non-meromorphic terms. 
According to our partial result, it seems that for a^^^ model with n > 3, the non- 
meromorphic terms do not add up to zero while for rf^^^ theory, they do cancel among 
themselves exactly and very nontrivially. 

There are many things for future works. To mention a few, the first thing 
is to generalize this method systematically to all loop order, which seems rather 
straightforward though the actual evaluation may not be easy. The second thing is 
to consider the multi-point correlation functions. The third thing is to accommo- 
date non-trivial boundary potentials maintaining the integrability of the model. Of 
course, the mysterious factor half and the non-vanishing of the non-meromorphic 
terms for al^^ model with n > 3 need much attention. 




[1] 



(40) 



'n 
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